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监测 模式 (跟踪 、 扫 
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W, T UAE BS OT FE 2k 


(天 气 和 设备 ， 


测 网 的 整体 效能 。 


本 文 认 为 ， 监 测 网 的 主要 任务 是 采集 足以 进行 碎片 常规 
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给 定 一 个 日 
任务 调度 流程 如 图 


拟 对 M 个 空间 碎片 进行 观测 。 首 
观测 弧 段 信息 其 次 ， 


弧 段 的 时 空 分 布 和 


监测 网 任务 调度 优化 模型 的 建立 
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设备 的 位 置 和 空间 碎片 的 轨道 信息 ， 计 算 可 
赋予 该 弧 段 的 监测 收益 ; 


最 后 ， 基 于 可 测 


监测 收益 ， 解 算 监 测 任 务 优化 问题 ， 并 输出 各 设备 的 观测 任务 表 。 
/ 标 信息 / /监测 网 信息 / 
l — 任务 规划 优化 解 算 
/ TAS I 
各 装备 观测 任务 列表 编排 


每 个 可 观测 弧 段 收益 


图 1 
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站 监测 网 信息 后 的 时 间 分 布 情况 。 
对 一 个 碎片 进行 
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占 在 同一 时 段 有 
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表 1 列 出 碎片 - 台 站 的 


其 中 时 间 单 位 r 


Ke 


段 中 可 以 监测 到 目 


标 即 可 记录 目标 
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信息 。 


假设 每 个 台 站 对 碎片 监测 了 时 间 长 度 为 的 子 弧 段 即 
始 监测 下 一 个 碎片 。 对 任务 规划 时 间 进 行 


可 测 碎片 信息 ， 可 


以 获得 如 表 2 的 收益 列 寻 


个 有 效 监 测 弧 段 (或 需要 的 最 短 监测 六 


固定 时 长 7 的 栅 格 化 处 理 ， 相 
RRR， 按照 测 站 编 


可 观测 弧 段 时 间 及 对 应 的 监测 收益 ， 
段 ) 的 时 间 长 度 ， 如 果 在 7 时 间 


完成 有 效 监 测 ， 此 后 ， 便 立即 开 
民 据 栅 格 化 的 时 间 和 
号 作业 者 序列 ， 作 


号 大 小 往 后 编 
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Ra ER IAF 弧 段 结 


za 
编号 ”编号 ” 始 时 间 ” 束 时 间 ne 
1 I 0.97 2.27 1-I 
1 II 1.87 3.47 LI 
1 IV Or 1.37 1-IV 
1 IV 2.57 3.97 1-IV 
2 I 2.57 4T 2-I 
2 II 1.37 2.17 2- I 
2 M 0.97 3.87 2- 
测 站 1 测 站 2 
2 2 个 台 站 对 4 个 碎片 的 可 测 弧 段 时 间 分 布 
业者 的 概念 将 会 在 2.2 节 中 正式 提出 。 表 中 有 注释 处 即 为 有 观测 弧 段 记录 。 
表 2 测 站 对 碎片 的 栅 格 化 收益 列表 
作业 者 
Ht 1 2 3 4 5 6 7 8 
I 1-I 2-1 
I LI 
M 2-II 2- 
IV 1-IV 
2.2 ”组 网 调度 优化 数学 模型 
基于 上 述 场景 ， 监 测 网 的 空间 碎片 任务 规划 优化 可 用 如 下 方程 表示 : 
N M 
MaxZ = XOY byti 5 (1) 
i=1 j=l 
其 中 ， 
S | 1, 作业 者 i 监测 到 碎片 ; 匹配 (可 观测 ) 
ij = 0, 
M 
Sas, sav ’ (3) 
j=l 
N 
Yo tes 12.8; ,M ) (4) 
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其 中 ，i 是 栅 格 化 后 的 监测 任务 序号 ， 
栅 格 化 时 间 段 个 数 的 乘积 
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测 站 地 理 位 置 、 设 备 观 测 时 段 空 


E 同 决定 ， 既 可 以 反映 碎片 的 
方 (用 户 ) 根据 碎片 近期 监测 弧 段 的 分 布 和 数据 质量 、 
s 闲 情况 等 


称 之 为 作业 者 ， 
; ERIAREN M; bi; 是 作业 者 i 监测 碎片 7 时 的 
要 性 ， 也 可 以 反映 测 站 的 重要 性 。 
所 需 轨 道 精度 、 距 离 上 次 观测 时 长 和 
素 综合 确定 。 当 监测 对 象 为 重点 关注 的 目 


El 
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j 户 可 能 要 求 只 要 有 机 会 
以 确保 该 目标 可 见 时 一 定 会 被 某 测 站 所 


监测 。 


就 对 该 目标 进行 跟踪 观测 ， 
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作业 者 的 总 数 为 N，N 


bij 可 由 
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监测 收益 设 为 一 


个 很 大 的 值 ， 


能 匹配 一 个 碎片 进行 
维 
轨 精 度 要 求 确定 ) 得 到 一 次 监测 ， 
等 情况 。 在 式 (3) 和 


监测 ， 


es 


首先 
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成 。 另 一 方面 ， 如 果 一 个 碎片 可 同时 被 多 台 设 备 监测 ， 


备 空闲 (无 指派 任务 )， 
A (1) 一 (4) FERH 
是 运筹、 优化 领域 研究 的 重点 ， 其 
前 ， 计 算 机 常用 的 成 熟 算法 
algorithm)" il Bertsekas 和 Orlin # 

展 而 来 的 各 种 优化 算法 。 其 中 ，LAPJV 


这 些 


台 对 偶 算 法 (primary dual algorithm, PD) 基础 上 


= Bae HT Bi fr RRUA 
指派 问题 ” (linear assignment problem， 的 数学 表达 。 
解法 主要 有 原始 对 人 
包括 Kuhn 和 Munkres @ 


等 人 的 拍卖 算法 


算法 


bi; 的 具体 确 
设备 监测 任务 执行 效能 等 进行 调整 ， 是 需要 更 进一步 研究 的 问题 。 
式 (4) 表示 一 个 碎片 在 该 段 时 间 
护 / 扩 展 为 主要 任务 的 监测 网 任务 规划 ， 
而 何 时 
(4) 的 约束 下 ， 将 确保 监 六 


遇 算 法 、 


定 方式 将 随 碎片 编目 性 
式 (3) 表示 一 个 作业 者 只 
内 最 多 被 监测 一 次 。 


能 表现 和 


以 碎片 编目 


收益 
个 碎片 的 监测 任务 完 


需要 确保 必要 的 碎片 在 一 定时 间 内 (根据 定 
由 何 设备 监测 则 须 综合 考虑 可 见 性 、 
任务 总 体 收益 最 大 且 


(紧迫 性 ) 


除 被 指定 的 设备 监测 外 ， 如 果 其 他 设 


行 监测 ， 


(auction algorithm)’ 


path algorithm, SAPA) 提出 的 一 种 高 效 的 求解 


2.3 LAPJV 算 法 的 改进 


法 。 


LAPJV 算法 在 相同 的 问题 规模 下 


顾 多 个 优 


1) 开辟 内 存 小 ， 


( 

(2) XIF IE E5 AS E E PE BD eT A EB H 
(3) 对 于 b;; 组 成 的 成 本 (或 收益 ) FRE B H 
但 大 规模 组 网 的 现实 场景 中 用 此 算法 
( 
(2 
(3 


1) 作业 者 总 数 NN 和 碎片 总 数 M 往 


) 数据 繁多 ， liege aia. 
度 非 常 


这 些 难 


) JE EFR D 
考虑 到 


a 
(1) 为 了 获得 成 本 和 矩阵 的 方 阵 表达 ， 引 入 虚拟 作业 者 或 者 虚拟 碎片 对 入 


往 不 相同 ; 


A, afi: 
代码 量 适中 ， 较 少 的 解 算 次 数 却 


有 较 高 的 解 算 速度 ; 


4 取信 敏感 度 低 。 
会 存在 如 下 三 个 问题 : 


法 进行 了 如 下 修改 : 


单纯 形 法 、 对 介 
等 人 的 匈牙利 算法 


[14-16 


本 


并 不 会 影 ! 


向 总 体 规 划 。 
LAP 一 直 


y, [10] 
法 等 。 


(Hungarian 


以 及 由 此 发 


(LAP solved by Jonker and Volgenant) 是 在 原 
结合 最 短 增 广 路 径 算法 


(shortest augment 


E 阵 进行 扩充 ， 表 
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达 式 如 下 : 
L L 
Max Z = XOX byty = BrxL: XLxL = 
fat 
bi © biL Til t TIL 
(L = Max{M,N}) . (5) 
br1 © brL| Zr © TLL 


(2) OAS HE BERN AT EENE T FERE SRE AE. PA ES 
元 素 和 零 元 素 可 以 通过 建立 BLxz 的 一 种 新 映射 进行 有 效 存储 : 首先 ， 删 去 无 效 的 零 行 和 
零 列 ， 包 括 因为 式 (5) 增加 的 部 分 ， 重 新 组 织 成 本 矩阵 减少 稀 朴 度 和 从 形式 上 缩小 矩阵 维 
ee a T 
REA BUS MRR ENAME Ba, AA EMM fe E e REKI LAPJV 
BATE A Ji GAS EY 4 aed JB BR VE. MO J A LAPJV pine LAPJV, 
ILAPJV) SHS BE 76 BOK AY SE Ta BE a a it FE RE BERI A ETE fi Pe OT o 


3 监测 网 任务 调度 优化 模型 的 改进 


3.1 ”线性 模型 的 不 足 与 改进 

在 线性 指派 模型 中 ， 无 论 下 一 个 碎片 的 位 置 如 何 ， 算 法 都 会 选择 使 全 局 收益 最 大 的 可 
测 碎片 ， 而 不 考虑 望远镜 从 当前 位 置 移动 到 监测 该 碎片 所 需 的 成 本 。 如 图 3 所 示 ， 碎 片 1 
(Debris 1) 被 监测 完成 后 ， 如 果 监 测 碎片 2 (Debris 2) 与 3 (Debris 3) 的 收益 相等 ， 则 可 能 
选择 望远镜 移动 成 本 更 低 的 碎片 2， 如 果 监 测 碎片 3 的 收益 更 高 ， 同 时 移动 望远镜 到 监测 碎 
Fr 3 的 时 间 成 本 也 更 高 ， 此 时 ， 应 综合 考虑 监测 收益 和 望远镜 移动 的 时 间 成 本 。 实 际 上 ， 当 
前 的 每 一 步 选择 都 会 影响 之 后 的 选择 ， 局 部 与 全 局 之 间 存 在 矛盾 。 

用 Ci; 表示 将 望远镜 (作业 者 i) 移动 到 监测 碎片 7 位 置 的 成 本 ， 而 这 一 成 本 取决 于 当前 
望远镜 与 待 测 碎片 的 位 置 ， 用 dil, k) 表示 碎片 无 ( 刚 完 成 监测 的 碎片 ) 与 碎片 7 之 间 的 角度 
视 距 ， 望 远 镜 移动 的 成 本 由 下 式 给 出 ; 


Cy =e (Yon 1b ° di(j, k I)» ) (6) 


i=2 j=l 


AP, c 是 移动 望远镜 的 单位 成 本 (望远镜 移动 0.1 弧度 所 需 的 时 间 代 价 )。 当 前 伺服 能 力 下 ， 
地 基 望 远 镜 远 距 离 移 动 与 近 距 离 移动 的 时 间 成 本 并 非 是 线性 增长 而 是 相对 接近 的 ; 但是， 为 
了 保持 卫星 姿态 的 稳定 ， 天 基 望 远 镜 的 移动 时 间 基 本 上 与 移动 距离 呈 线 性 关系 。 无 论 是 天 基 
还 是 地 基 的 跟踪 监测 ， 望 远 镜 移动 成 本 均 是 一 个 值得 考虑 的 因素 。 望 远 镜 移动 成 本 的 具体 确 
定 方式 及 其 对 监测 任务 规划 的 影响 ， 是 下 一 阶段 需要 研究 的 问题 。 一 般 而 言 ， 描 述 成 本 的 数 
值 需 与 收益 值 相 匹配 。 
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测 站 


À 工作 者 序列 


碎片 1: 
结束 跟踪 


d(j,k) 


六 片 1 结束 时 
卒 片 2 的 位 置 


碎片 1: 
开始 跟踪 


碎片 1 结束 时 
碎片 3 的 位 置 


图 3 考虑 移动 成 本 调度 示意 图 


结合 式 (6) 与 (1)， 综 合 考 虑 监测 收益 和 望远镜 移动 成 本 ,任务 调度 优化 目标 函数 就 


变 成 
Max Z’ = Benefit — Cost = 
N M 
> >， C D 1,kdi( (j,k j) a (Hi = 1 时 ， C 一 0) (7) 
站 二 


其 余 的 约束 条 件 均 与 式 (2) 一 (4) 相同 。 此 非 线性 问题 与 Dantzig 等 人 提出 的 旅行 商 问 
题 (travelling salesman problem, TSP)” 有 相似 之 处 。TSP 一 般 描 述 如 下 : 为 推销 员 找 到 
从 给 定 城市 出 发 ， 访 问 指定 城市 组 中 的 每 个 城市 ， 然 后 返回 原始 出 发 点 的 最 短路 线 (旅程 )。 
3.2 ” 非 线性 模型 与 TSP 对 比 

TSP 的 最 优 解 问题 已 经 被 证 明 是 一 个 NP-Hard 问题 ”， 其 时 间 复 杂 度 O(n!) 为 指数 型 
时 间 复 杂 度 。 目 前 解决 TSP 主要 有 两 个 方向 :一 个 是 耗 时 长 但 精确 的 最 优化 算法 ， 另 一 个 
ese, HEN: HEME. HEN he 
等 。 如 果 能 够 接受 微小 偏差 ， 启 发 式 算法 会 在 较 短 时 间 内 获得 较 好 的 局 部 最 优 解 ， 且 解 的 结 
果 十 分 接近 全 局 最 优 解 。 

NB a te er i 

(1) 碎片 与 测 站 性 质 不 同 ， 无 法 直接 建立 一 个 距离 方 阵 ， 可 能 需要 扩充 成 极 大 维度 的 方 
阵 ， hela glans a Seca ope niente 

(2) 不 是 所 有 的 测 站 都 有 可 观测 的 碎片 ， 且 需要 检验 碎片 是 否 已 被 唯一 观测 ; 

(3) 解 算 路 径 不 是 闭环 。 

考虑 到 这 些 问 题 并 经 过 小 规模 的 测试 ， 笔 者 认为 仿生 类 智能 算法 并 非 适 宜 的 算法 ， 而 
k-opt 算法 能 较 好 地 解决 TSP 问题 。 
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城市 
碎片 
o 
Zi 
Jasta 
作业 者 
E m 


4 ”考虑 移动 成 本 非 线 性 模型 和 TSP 规划 示意 图 


3.3 k-opt (k =2、3) 算法 的 改进 


准 
更 
大 


分 


k-opt 算法 ， 也 称 k-exchange 算法 ， 它 的 策略 是 在 顺序 交换 准则 、 可 行 性 准则 、 正 收益 
则 和 不 相交 准则 作用 下 ， 在 随机 初始 回路 的 基础 上 ， 通 过 不 断 地 更 换 K 个 链 路 连接 进行 
新 ， 直 到 其 不 能 产生 更 短 的 路 径 长 度 。 大 的 取 值 越 大 ， 最 后 所 得 解 是 最 优 解 的 可 能 性 越 
[21,22] 

然而 ，k-opt 算法 的 时 间 复 杂 度 是 O(n*)。 研 究 表 明 2-opt 算法 搜索 并 解决 TSP 问题 十 
9 效 “” ， 文 章 将 在 实验 部 分 探究 2-opt 与 3-opt 算法 实验 结果 的 差异 性 。 为 了 使 k-opt 更 


适合 非 线性 模型 的 解 算 ， 对 它 做 出 四 处 修改 。 


在 此 基础 上 减 去 望远镜 移动 成 本 以 生成 初始 解 。 


大 


1% 


4 


(1) 根据 k-exchange 算法 主要 思想 进行 局 部 交换 ， 但 不 严格 遵循 其 解 算 流 程 。 
(2) 采用 贪心 算法 (greedy algorithm, GA) 给 作业 者 选择 收益 最 大 且 不 冲突 的 碎片 ， 并 


(3) (EFA Pit REM id ay slid ae A OUI fa, SET] ASE eee ar A I Be ez EC 
小 排序 再 搜索 。 
(4) 尝试 交换 进 新 的 碎片 时 ， 会 使 用 特殊 的 检查 函数 加 快 收敛 速度 ， 保 证 作业 者 与 碎片 


匹配 唯一 性 ， 帮 助 我 们 只 想 尝试 最 有 可 能 增加 效益 的 作业 者 -碎片 组 合 。 例 如 一 个 较 高 收益 


片 已 被 尝试 并 失败 时 ， 另 一 个 低 收益 碎 片 则 不 会 尝试 。 


实验 结果 


4.1 ILAPJV 与 k-opt (k =2、3) 算法 性 能 比较 


F, 


为 了 测试 在 数据 量 较 多 情形 下 ILAPJV 算法 的 性 能 与 kopt (k =2 3) 算法 的 性 能 
设计 了 两 个 任务 调度 优化 实验 。 第 一 个 实验 中 ， 给 定 4 个 测 站 、1000 个 空间 碎片 的 轨 


道 参数 ，3 d 内 的 可 监测 弧 段 达到 46000 个 之 多 ， 每 个 弧 段 的 监测 收益 为 2~32 的 数 ， 第 二 
个 实验 中 ， 给 定 4 个 测 站 、5000 个 空间 碎片 的 轨道 参数 ，3 d 内 的 可 监测 弧 段 达到 220000 


之 


多 ， 每 个 弧 段 的 监测 收益 为 0~800 的 数 。 设 有 效 监测 弧 段 的 长 度 为 4 mins WEA A LE 


分 析 中 ， 某 些 设备 对 某 个 目标 单 次 最 长 可 连续 监测 50 min， 从 轨道 编目 角度 而 言 ， 只 须 监 


测 


4 min; 但 如 果 该 目标 为 重点 目标 ， 则 可 通过 设置 极 高 收益 对 其 持续 监测 。 实 验 环境 为 


4 期 


钱 阳 波 ， 


Intel(R) Core(TM) i5-10400F CPU @ 2.90GHz. 


分 别 进行 120、240 和 480 min 的 监测 
其 中 精度 指标 为 革 介 


时 间 、 总 收益 和 精度 ， 


作业 规划 ， 


算法 的 性 能 指标 为 任务 调度 优化 解 算 
法 规划 的 任务 收益 与 ILAPJV 算法 规划 的 任务 
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收益 之 比 ， 以 百分比 表示 。 第 一 个 实验 中 ， 将 ILAPJV 算法 、2-opt 算法 与 3-opt 算法 解 算 


线性 模型 的 性 能 指标 进行 对 比 ， 见 表 3; 表 4 是 第 


表 3 实验 1 优化 算法 性 能 指标 比较 


二 个 实验 的 结果 。 


算法 解 算 时 间 /s 收益 精度 /% 
120 240 480 120 240 480 120 240 480 
ILAPJV | 0.006 0.029 0.133 | 256.02 608.53 1452.87 100 100 100 
2-opt 0.017 0.044 0.120 | 256.02 608.23 1443.86 100 99.951 99.380 
3-opt 4.214 8.383 52.666 | 256.02 608.43 1451.15 100 99.984 99.882 
GA 0.005 0.007 0.011 | 256.01 601.42 1392.53 | 99.996 98.832 95.847 
TE: 120. 240. 480 分 别 表示 规划 监测 时 长 ， 单 位 为 min。 


表 4 实验 2 优化 算法 性 能 指标 比较 


算法 解 算 时 间 /s 收益 ta SE /% 
120 240 480 120 240 480 120 240 480 
ILAPJV | 0.053 0.128 0.364 | 12098.02 18794.26 30670.76 100 100 100 
2-opt 0.024 0.049 0.133 | 12096.44 18780.28 30656.60 | 99.987 99.926 99.954 
3-opt 0.288 5.085 46.083 | 12097.38 18793.46 30666.95 | 99.995 99.996 99.988 
GA 0.021 0.021 0.029 | 11955.50 18718.34 30143.15 | 98.822 99.596 98.280 


YE: 120. 240. 480 分 别 表示 规划 监测 时 长 ， 单 位 为 min. 


由 表 3 和 表 4 可 知 : 


(1) 两 个 实验 中 ILAPJV 


率 ， 且 收益 值 bij 的 大 小 对 算法 解 算 时 间 影 响 小 。 


(2) 2-opt 算法 和 3-opt 算法 都 有 不 俗 的 表现 ， 两 个 实验 的 解 算 精 度 者 


TN 


法 性 能 表现 优良 ， 对 于 较 大 数据 量 的 问题 有 很 好 的 解 算 效 


在 99% 以 上 。 实 


验 2 中 采用 2-opt 算法 解 算 优化 问题 所 需 的 时 间 比 ILAPJV 算法 更 少 ， 且 解 算 精 度 接近 或 


È F 99.9%。 


其 解 算 时 间 代 价 将 难 


采用 3-opt 


的 百倍 之 多 。 


法 时 ， 其 精度 比 


2-opt 算法 略 有 提高 ， 但 是 随 着 问题 


(3) GA 的 解 


综合 以 上 三 点 可 知 ， 采 月 


效率 最 高 ， 但 解 人 


确 ， 但 其 解 算 时 间 成 本 较 高 ， 不 适 


精度 不 其 理想。 
H ILAPIV 能 高 效 地 解决 线性 优化 问 
于 实时 任务 规划 优化 。 


以 忽略 ，240 min 的 任务 规划 所 需 解 算 时 间 是 对 


2-opt 算法 解 算 涉及 望远镜 移动 成 本 的 非 线 性 优化 问题 。 


以 实验 2 中 120 min 任务 规划 为 例 ， 给 定 4 个 测 站 、5 000 个 空 


位 成 本 c= 10、50 时 非 线 性 模型 可 监测 120 个 碎片 ，c = 100 可 监测 117 个 碎片 ， 表 明 移 


规模 变 大 ， 


< 用 2-opt 算法 


题 ， 尽 管 3-opt 结果 更 精 
因此 ， 采 用 兼顾 精度 和 速度 的 


所 需 TiTi 时 间 


x 间 碎片 ， 线 性 模型 和 单 
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动 成 本 会 影响 碎片 监测 数量 。 考 虑 每 个 监测 弧 段 长 度 为 4 min， 每 个 测 站 在 120 min 内 最 多 
可 监测 30 个 碎片 ，4 个 测 站 最 多 可 监测 120 个 碎片 ， 因 此 ， 监 测 设备 得 到 充分 利用 。 并 且 ， 
规划 的 监测 收益 最 大 。 如 图 5 所 示 ， 图 中 c = 0 表示 用 ILAPIV 解 算 线性 模型 后 的 监测 碎 
片 总 数 ，c 关 0 表示 用 2-opt 解 算 非 线性 模型 后 的 监测 碎片 总 数 。 


国 四 <-0 BD —-1 Ml Bll] 10 
m20 M210 Mp} c50 [2] 100 
500 


400 
300 
200 
120 
120 103 | | 

1 l jii 
, 

120 


监 ; ae /min 


碎片 总 数 


监测 


© 
© 


YE: [1] 表示 实验 1，[2] 表示 实验 2. 
图 5 监测 碎片 总 数 随 不 同 规划 模型 和 时 长 变化 图 


4.2 ”组 网 监测 站 布设 和 空间 碎片 选取 

组 成 空间 碎片 监测 网 时 ， 我 们 选用 7170 个 NORAD 编目 库 中 的 空间 碎片 ， 其 ed 
道 ” 的 产生 采用 陈 俊 宇 等 人 提出 的 基于 TLE 数据 的 编目 轨道 改进 方法 ， 通 过 武汉 大 学 空 
间 态 势 感知 团队 开发 的 效能 软件 设置 地 面 测 站 并 进行 碎片 可 见 性 分 析 。 监 测 网 由 遍布 全 球 
的 200 个 地 面 测 站 组 成 ， 测 站 位 置 选 自 国 际 GNSS 服务 (International GNSS Service, IGS) 
发 布 的 SNX 文件 ， 如 图 6 所 示 。 


设 定 每 个 测 站 均 配 置地 基 跟 踪 式 望远镜 


测 站 设备 参 3 

ao o aa ERB MILA 5. 
探测 敏感 度 20 mag 7 显示 了 效能 软件 生成 仿真 观测 弧 段 的 
方位 角 0~360° 方法 : 在 输入 目标 的 精密 轨道 、 测 站 的 站 址 和 
高 度 角 0~90° 设备 信息 、 仿 真 起 止 时 间 等 基本 参数 后 ， 根 据 
有 效 监 测 弧 段 长 度 时 间 60 s 设 定 的 步 长 对 目标 轨道 进行 内 插 ， 在 每 一 个 内 
ee 2n 插 时 刻 计算 目标 与 测 站 在 惯性 系 中 的 位 置 ， 根 


据 观 测 站 观测 范围 ， 计 算 目 标的 可 见 性 ， 若 目 
标 相对 于 测 站 可 见 ， 则 计算 目标 观测 值 ， 生 成 测 轨 弧 段 。 
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图 6 200 个 地 面 观测 站 组 成 的 空间 碎片 监测 网 


目标 与 测 站 位 置 计算 


可 见 性 计算 


观测 弧 段 生成 


望远镜 观测 弧 段 


图 7 效能 软件 对 可 见 弧 段 的 仿真 
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4.3 ”全 球 组 网 监测 仿真 结果 


根据 
te 


足 最 


利用 以 上 地 面 网 监测 碎片 ， 监 测 作业 时 间 为 2021 年 3 月 


41 卷 


7 日 至 2021 年 3 月 9 日 ， 


个 碎片 的 监测 收益 设 定 为 0~100 之 间 的 随机 数 。 设 有 效 监 测 弧 段 的 时 间 长 度 为 60 s， 这 是 
编目 目标 轨道 更 新 和 新 目标 初 轨 确 定 / 弧 段 关 联 等 算法 对 弧 长 要 求 而 设置 的 。 进 行 4h 
测 任务 优化 ( 方 阵 为 Bag000x4g000)， 同 一 设备 对 同一 目标 单 次 最 长 连续 监测 180 s， 满 
长 观测 时 长 后 即使 可 见 也 不 再 观测 。 计 算 所 需 时 间 和 监测 收益 如 表 6 所 示 ， 结 果 显 示 : 

表 6 4h 规划 任务 的 仿真 结 
ap Ky 景 WAY 则 A 
| ah Ye) K/s 总 收益 人 ase o 
ILAPV, c=0 12.051  289399.07 40.88 
2-opt, c=1 162.071 285 333.79 40.68 
Qopt, c=10 164.026 257953.53 40.70 
Qopt, c=50 160.941 147811.66 35.73 


效率 


(1) 使 用 ILAPJV 算法 解 算 线 性 问题 (不 考虑 望远镜 移动 成 本 ) 的 时 间 为 12.051 s， 解 算 


很 高 ; 


(2) 设 望远镜 移动 成 本 为 每 弧度 c = 10， 使 有 


可 以 控制 3 min LAA; 
(3) 时 长 48 的 观测 方案 
共有 2931 均 占 全 部 碎片 40% 以 上 。 


模型 


总 的 设备 利用 率 。 H 
进行 任务 规划 ， 


设 为 


TE Hr Ay SL 


3 min, 


E 分 析 表 明 ，200 个 台 
HAS EY TA] 1 min 时 ， 设 备 利 
All FA EZ t 


， 利 用 非 线性 模型 共有 2918 个 碎片 被 观测 了 


台 站 中 只 


E 和 线性 模型 


14 158 个 台 


站 有 可 见 碎片 可 以 观测 ， 


此 时 ， 设 备 的 时 间 利 用 率 基本 上 是 1 min 时 的 3 倍 。 


7000 多 碎片 ， 监 测 资 
求 进 行 任务 规划 。 
目标 的 监测 任务 规划 。 


主要 


行 监 
于 碎 
Te Ar 


监测 


有 的 碎片 有 望 被 观测 一 次 ， 因 


测 作业 ， 显 然 可 提 
片 数 量 比较 紧缺 时 ， 设 备 利 
作业 任务 完成 情况 发 生变 化 时 ， 可 及 时 


结合 这 两 个 因 


高 设备 利 


j 率 看 起 来 不 高 。 


H 2-opt 解决 非 线 性 问题 的 时 间 为 164.026 s, 


一 次 ， 而 线性 


从 而 影响 了 


为 此 ， 我 们 将 栅 格 时 间 
仍 可 分 别 观测 2918 个 和 2931 个 碎片 ， 


考虑 到 该 实验 中 有 200 个 台 站 而 只 有 


供 源 比较 丰富 ， 再 考虑 到 我 们 按 每 个 碎片 只 需 观 测 一 次 这 一 
素 ， 可 以 认为 ， 本 文 提出 的 


法 较 好 1 


男 一 方面 ， 在 监测 资源 比较 丰富 时 ， 


1, 2 


以 设 站 


地 完成 了 以 编目 维 
n 
且 不 会 影响 总 体 规划 ， 可 


基本 监测 需 
护 为 
自行 或 被 指派 进 
， 当 设备 资源 相对 


j 率 会 较 高 。 


实际 上 ， 


40% 碎片 这 一 结果 ， 如 果 接 下 来 短 时 间 内 不 需要 


4.4 其 他 任务 说 明 


型 、 考 虑 移动 成 本 的 非 线 性 


总 结 模型 


局 最 


优 场景 。 


对 于 异常 目 


标的 快速 监测 和 习 


进行 新 的 规划 。 
次 对 这 些 碎片 j 
此 ， 可 以 较 好 地 完成 编目 库 维 
上 述 算法 测试 结果 表明 ，2-opt 算法 可 用 于 大 规模 监测 


要 目标 的 连续 监测 ， 


PACS 
最 后 ， 


护 的 目标 。 


4 和 算法 ， 笔 者 将 上 述 工作 转化 成 一 个 碎片 监测 任务 规划 软件 。 
模型 的 测试 ， 还 是 大 规模 场景 的 仿真 ， 都 8 


网 的 近 实 时 监测 


也 认 使 用 常规 观测 的 全 


法 的 近 实 时 规划 能 力 ， 当 


根据 4 h 内 观测 了 


行 监 测 ， 在 10h 内 ， 所 


任务 规划 优化 。 


无 论 是 线性 模 


可 通 


过 图 8 所 示 软 件 右 侧 
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可 视 化 界面 的 交互 功能 完成 
| E - 
Input fle [FA Browse. 5o Debris 10 [IED Station Number [1 
Bead input fle ea 
raona [ -3 
SUse manual benefit for fast moving objects 115 Set Time |2452214.510995 Set Elevation [0.261440 SetAzimuth [6.246720 $ 
Ciena tn 7 ton [ROHS 4 Cx] 
r hasan eS 121 厂 Track Entire Pass Cancel 
e i 
E em | || 1 
图 8 ”规划 软件 操作 界面 
为 了 使 某 一 异常 碎片 被 纳入 监测 列表 ， 可 在 右 侧 甘 特 图 中 使 用 鼠标 右 击 打开 功能 窗 
选择 “强制 监测 ” 如 果 想 使 重要 目标 实行 连续 监测 ， 可 在 功能 窗口 选择 “跟踪 全 程 ”。 前 
者 其 实 是 通过 强行 赋值 碎片 高 收益 值 使 其 不 得 不 被 监测 实现 ， 后 者 通过 强制 选择 可 见 弧 段 
实现 。 此 外 ， 其 过 境 时 间 可 能 不 足 一 个 最 短 有 效 观 测 时 段 ， 如 图 9 中 的 碎片 A， 或 跨越 最 
短 有 效 时 段 ， 如 雁 片 B。 如 果 需 要 监测 这 类 碎片 ， 会 丢失 一 小 部 分 网 格 时 间 ， 故 建议 在 效益 
极 高 时 调整 网 格 。 实 际 上 ， 大 碎片 能 被 设备 监测 ， 在 它 的 (其 ) 短 弧 初 轨 确 定 、 编 目 工 作 中 ， 
并 非 需 要 对 完整 弧 段 进行 监测 ， 从 这 一 角度 而 言 ， 栅 格 化 损耗 不 会 十 分 普遍 ， 少 量 的 栅 格 化 
损耗 可 以 被 接受 。 此 处 栅 格 化 损耗 指 : 本 来 纳入 监测 工作 可 完整 监测 的 一 段 有 效 时 间 ， 由 于 
特殊 需要 并 没有 完全 利用 ， 如 图 9 栖 色 部 分 表明 的 空闲 时 间 。 
— 碎片 B 
A 
碎片 i | BEEE 
HAH | 一 | | | 
1 I | 作业 者 序列 
VE: 蓝 色 横 条 为 碎片 可 监测 时 段 ， 黄 色 横 条 是 涉及 的 有 效 观 测 时 段 ， 橙 色 部 分 是 损耗 时 间 。 
9 ”特殊 碎片 规划 任务 说 明 
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5 结束 语 


针对 多 任务 场景 下 空间 碎片 组 网 监测 任务 的 调度 优化 问题 ， 我 们 分 别 建立 了 不 考虑 望 
远 镜 移动 成 本 的 线性 模型 和 考虑 望远镜 移动 成 本 的 非 线性 模型 ， 后 者 更 真实 地 描述 了 实际 
监测 问题 。 在 分 析 线 性 问题 特性 的 基础 上 ， 改 进 了 LAPJV 算法 ， 实 现 了 快速 精准 的 优化 问 
题 全 局 最 优 解 。 对 于 非 线性 模型 ， 利 用 改进 的 2-opt 算法 ， 可 实现 局 部 最 优 解 的 高 效 解 算 ， 
解 算 结果 与 全 局 最 优 解 具有 较 好 的 一 致 性 。 这 两 个 模型 具有 非常 现实 的 意义 。 算 法 性 能 测 
试 结果 表明 ， 改 进 的 两 种 算法 可 用 于 解决 监测 网 任务 调度 时 数据 量 繁多 、 实 际 有 效 观 测 值 稀 
少 等 问题 ， 上 共有 较 高 的 解 算 精度 和 效率 ， 可 用 于 大 规模 监测 网 的 近 实 时 监测 任务 的 规划 优 
化 。 


虽然 在 短 时 间 内 ， 巨 量 的 空间 碎片 无 法 得 到 清除 ， 但 一 方面 ， 以 空间 天 气 监测 和 研究 为 
核心 的 监测 设备 正 逐 步 拓 展 到 空间 碎片 领域 ”。 我 国 天 地 一 体 化 空间 目标 监测 系统 正在 积 
极 建设 中 ， 全 球 范围 内 的 监测 网 络 日 益 庞大 ， 国 际 合作 趋势 显现 。 另 一 方面 ， 频 发 的 碎片 事 
件 和 太空 热点 活动 等 使 空间 监测 任务 更 趋 复杂 多 样 ， 对 时 效 性 提出 更 高 的 要 求 ， 监 测 网 整体 
效能 优化 成 为 重要 的 技术 手段 。 在 这 样 的 情形 下 ， 基 于 多 任务 场景 下 的 线性 和 非 线性 的 规划 
模型 和 解 算 算法 可 以 为 近 实 时 的 调度 要 求 提 供 一 种 有 效 的 解决 方案 。 


zi 
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A Scheduling Optimization Method for Space Debris 


Monitoring Network in Multiple Monitoring Task Scenarios 


QIAN Yang-bo!, XIA Sheng-fut, SONG Zheng-ji?, SANG Ji-zhang! 


(1. School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China; 2. Beijing Institute 
of Spacecraft System Engineering, Beijing 100094, China) 


Abstract: Orbit measurement data collected by a space debris monitoring network is the 
basis of debris orbit cataloguing. With the vast number of debris objects but limited moni- 
toring facilities, data collection strategy and task programming are the key technologies to 
bring into full play of monitoring facilities and improve cataloguing capability and accuracy. 
The monitoring tasks include routine monitoring, important object monitoring and monitor- 
ing of emergent events. For the multi-task scheduling optimization problem of a monitoring 
network, this paper proposes a linear assignment model and a nonlinear assignment model 
considering move cost, respectively, in which the monitoring benefit is defined as the ob- 
jective function, and uses the improved LAPJV algorithm and improved 2-opt algorithm to 
solve them. The optimization experiments with ground-based simulated debris monitoring 
networks are carried out, in which the improved LAPJV algorithm and 2-opt algorithm have 
a solution time of 12.051 s and 162.071 s, respectively, and the total benefit is 289 399.07 and 
285 333.79, respectively in term of the linear model and the nonlinear model with 4-hour task 
scheduling for 200 stations and 7170 debris. The final results show that both the algorithms 
are able to solve near real-timely the monitoring task programming, and can be used as an 


effective solution to the network task scheduling optimization. 


Key words: space debris; monitoring and maintenance; scheduling optimization/programming; 
LAPJV; k-opt 


